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Abstract. - We study the effect of fractal initial conditions in closed reactive systems in the 
cases of both mobile and immobile reactants. For the reaction A + A — ► A, in the absence 
of diffusion, the mean number of particles A is shown to decay exponentially to a steady 
state which depends on the details of the initial conditions. The nature of this dependence is 
demonstrated both analytically and numerically. In contrast, when diffusion is incorporated, 
it is shown that the mean number of particles (N(t)) decays asymptotically as t~ df ^ 2 , the 
memory of the initial conditions being now carried by the dynamical power law exponent. The 
latter is fully determined by the fractal dimension df of the initial conditions. 



Introduction. - In classical approaches to reactive processes mean-field (MF) theories 
have been useful to predict many non-trivial dynamical and steady state properties such as 
multistability, periodic concentration oscillations, chaotic motion etc. As long as fluctuations 
in concentration and occupation number space remain small, MF approaches provide a suitable 
description of the process. This is the case in many practical situations, where inhomogeneous 
fluctuations can be efficiently removed from the system by continuous external stirring of the 
reactants or sufficiently fast internal diffusion pQ. 

On the other hand, considerable interest has been devoted to reactive processes taking 
place on low-dimensional supports in recent years. In such systems, external stirring may 
prove difficult, and the internal diffusional mixing may not be sufficiently fast to compensate 
the effect of the restricted support geometry [2]. Correlated spatial fluctuations initially 
present in the system or induced by the interplay between the chemistry and the spatial 
characteristics of the support may then give rise to anomalous temporal behaviour ; 3: and 
even deviations from the MF steady state 0] . In particular, the outcome of the reaction may 
be strongly influenced by the dimensionality of the support. 

In recent years, the important role of dimensionality was recognized to be an issue common 
to many statistical processes as well |4I5| . Fractal supports, which can be produced to display 
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a continuum of dimensionalities, have been used widely in experimental |S1[7] and theoretical 
studies [S] to illustrate the new effects arising in the above context. A second factor of 
importance in the interplay between local dynamics and transport processes is the form of the 
initial conditions |2|. Depending on whether they are homogeneous, short-ranged or correlated 
on a long range they may favour quick mixing by diffusion, delays in such a mixing resulting 
in segregation or, in the limit of immobile reactants, nonergodic behaviour altogether. The 
objective of the present work is to analyze the evolution of a low dimensional system subject 
to initial conditions corresponding to a fractal spatial distribution of reactants in the lattice. 
In particular, the role of the fractal dimension df of the initial distribution will be assessed. 

The specific system we shall analyze here is the reaction A + A — > A |1UU11| . The interest 
in the dynamics of this reaction on fractal supports first arised in connection with exciton 
fusion experiments on fractal percolation clusters [Sl[7|- While previous work has treated the 
case of random homogeneous initial conditions |12| and the general inhomogeneous case |13| . 
we focus here on the specific case of an inhomogeneous fractal distribution. We consider the 
cases where a) the reactants A are immobile on a fractal subset of a ID lattice and b) the 
reactants A are initially placed on the fractal subset, but for t > they can diffuse and react 
throughout the entire ID lattice. We show, analytically and numerically, that for immobile 
reactants the steady state reactant concentration does not depend only on df but also on the 
details of the fractal subset. In contrast, for diffusing reactants we show that the memory of 
the initial fractal distribution is carried by the anomalous dynamical exponent which is fully 
determined by df. 

Dynamics on fractals: immobile reactants. - To investigate the effect of fractal initial 
conditions on reactive systems with immobile reactants we consider the reaction A+ A — » A on 
two different Cantor-like sets C\ and C^. Set C\ is obtained by iteration of the segment 1110, 
while set C2 is formed by iteration of the segment 1101. After the n-th iteration, the lattice 
size L defined for both sets contains 4™ sites. In both sets, the number of ones and zeros is 3" 
and 4™ — 3™ respectively. We consider the " 1" sites as the active sites (sites where reactions 
can take place) and the "0" sites as inactive, or empty sites. Both sets are deterministic 
fractals with df — log(3)/ log(4). For convenience, we shall number the lattice sites from 1 to 
L, say from left to right, regardless of whether they are "0" (inactive) or "1" (active) sites. 

Next, we fill the "1" sites of set C\ or set C2 with particles A and let the reactions proceed 
according to the following Monte Carlo (MC) algorithm: at each time step At, a lattice site 
i is randomly chosen. If i and a randomly chosen neighbour i ± 1 are occupied, the particle 
A at site i is removed from the lattice with probability k^. In our model, the site choice is 
unbiased, i.e. sites i + 1 and i — 1 are chosen with equal probability. On the other hand, fcfj is 
the conditional probability of reaction at each time step given that two neighbouring sites are 
occupied and plays thus the role of a local reaction rate. We choose fixed boundary conditions 
by introducing two additional sites and L + 1 at the boundaries and specifying them as "0" 
sites. This particular choice of boundary conditions is selected for convenience and does not 
a play any important role in the large L limit. Finally, we set At = 7; (*) 

We are interested in the time evolution of the mean particle number {N(t)} and the 
associated mean concentration 9(t) :— (N(t))/L. Clearly, with the above choice of initial 
conditions, (iV(0)) diverges while 6>(0) vanishes in the limit L -> 00 as L( lo s( 3 )/l°s(4))-i = 
L t~ r . In order to obtain a well posed problem, we shall therefore consider the case of a 
finite system, as it is the case in experimental situations involving, for instance, mesoscopic 
scale devices (micelles, single crystallographic faces of a solid catalyst, etc.). Our objective 
will be to see whether some generic trends will nevertheless show up for long times and/or 



( 1 )This choice implies that all lattice sites are scanned once on average after one time unit. 
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large sizes L. We expect that the final number of particles (iV(oo)) will strongly depend 
on the specific form of the initial condition due to the lack of mixing, since particle islands 
evolve independently from each other. Thus, the spatial correlations present in the initial 
distribution will propagate in time. 

The function (N(t)) can be obtained heuristically by means of the following observation 
|14j : a string of fc consecutive sites (fc- tuplet) can be destroyed either by internal reaction 
events between particles inside the tuplet or by reaction between the particles at each edge 
site of the tuplet and particles sitting at occupied nearest neighbour sites outside the fc— tuplet. 
The latter events require the existence of a k + 1-tuplet. The dynamics of the mean number 
of fc— tuplets M^\t) (averaged over an ensemble of identical lattices) in a fractal set of size 
L = 4™ is given by the following hierarchy of equations |12II13| : 

^Mi n \r) = -(fc - l)M< n >(r) - (r), fc = 1 • • • k max , (1) 

where k max is the size of the largest fc-tuplet and (r) = for fc > k max . The reaction 
rate kn has been absorbed in the adimensional time variable r = k^t. The general solution 
of eqs. JIJ depends strongly on the details of the fractal set. It reads 

k — k 

Mt\r) = e-^J2 {e ~ T ~ iy M<$ B (0). (2) 

s=0 



For sets C\ and C2, k max — 3. The mean number of particles (singlets) is then 

(N(r)) = M[ n \r) = £ ^(0) + (M<">(0) - M<"»(0)) c - + e" 2 ^. (3) 



s=0 

In particular, the initial tuplet distribution for set C\ is given by the number of particle islands 
(= 3™ _1 ) times the number of fc-tuplets contained by each island. In this case, there are three 
singlets, two doublets and one triplet in each island, i.e. (0) = 3 • 3™ _1 = 3™, (0) = 
2 • 3"' 1 and M^ n) (0) = 1 • 3 n_1 . Asymptotically, we have 

(WM) = I j(^M£(0) = f = ffM = 5.3»-. (4) 

This gives a survival factor 77(71.) := (N (00)) / (N (0)) = h, i.e. the number of particles drops 
to half the initial value regardless of the lattice size. As emphasized by eq. 10)), each island 
yields asymptotically a mean number of particles equal to |. In contrast, a simple-minded 
combinatorial counting giving each final state of the island (ASA),(SAS),(SSA),(ASS) the 
same statistical weight yields the wrong factor |. This reflects the nonergodicity of the 
system, implying that the number of statistical paths leading to each steady state is different. 

For set C2 , the initial distribution is slightly more complex and contains a variety of island 
sizes. We have M^ n) (0) = 3™, M 2 (n) (0) = |(3" - 1) and M| n) (0) = ±(3™-* - 1). Using 

again eq. Q, we obtain (N(oo)) = 3™ + | and 77(71) = j% + ^5— ■ In the large L limit, 
77(71) -> ^ w 0.583. The size distribution of the islands plays a crucial role to determine the 
number of surviving particles. Particles in smaller islands have a higher survival expectancy. 
Therefore, more particles survive in set C2 than in set C\. 

For set C%, the combinatorial argument is again based on the size distribution of the islands, 
which is connected with the tuplet distribution through the equation i]^\t) = M^\t) — 
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2M™(r) + mI%(t) 13. We then have l[ n \0) = 4" } (0) = K3"- 1 + 1) and 4' l) (0) = 
l^gn-i _ According to the combinatorial counting, islands of size 1 and 2 will yield one 
particle asymptotically, whereas island of size 3 yield | particles on average. This is again 
wrong since, as we know, each three-particle island is reduced to | particles on average. 

The time evolution of (JV(r)) is easily computed for both fractal sets by substituting the 
corresponding expressions for the fc-tuplet distributions into eq. ©. In the long time limit, 
the dominant term describing the decay to the steady state is proportional to e~ T (recall that 
t = kat). Thus, the information on the initial distribution is contained in the coefficient 
of the dominant term rather than in the relaxation time fc^ 1 given by the exponent of the 
dominant term. In particular, this means that fractals with different df may relax at the same 
speed into the steady state. We therefore conclude that in the immobile reactant case neither 
the dynamics nor the steady state are suitably characterized by df. 

To confirm our results, we have performed MC simulations over 2 • 10 3 statistical runs on 
a lattice with L — 4 5 = 1024 sites. For both fractal sets, the asymptotic concentrations 6(oo) 
and the dynamics dictated by eq. © agree very well with the simulations (fig. 0) . 

Let us compare the previous results for inhomogeneous initial conditions with the case of 
a lattice containing only "1" sites. A lattice of length L = 4™ can be regarded as an iteration 
of, say, the segment 1111. The initial fc-tuplet distribution is given by M<; n) (0) = 4 n - fc + 1 
with k = 1, ..,4". In the large L limit, this yields (iV(oo)) = E^Lo^ 4 ™ ~ s )/ s! ~ e" 1 (4" - 1) 
and 77(71) w e _1 (l — 4~ n ), i.e., 77(71) — > er 1 « 0.367 as L~ x . As expected, the homogeneous 
system is characterized by a lower survival factor than any Cantor-like set for all values of L. 

A comment on the new features brought in by the fractal initial conditions analyzed above 
is now in order. The initial inhomogcneities imposed by sets C% or C2 decouple the dynamics 
of different parts of the system and decompose it into smaller, homogeneous subsystems. In 
principle, they have the same effect as the reaction-induced inhomogeneities, i.e. lowering the 
number of active sites. However, the interest of considering inhomogeneous initial conditions 
lies in the fact that, while the initial spatial correlations range over the whole system size, 
chemically-induced correlations are short-ranged; the latter may only develop between sites 
initially belonging to the same island. Due to the absence of diffusion or any other randomizing 
mechanism, a detailed memory of the initial spatial structure is carried by 6(t) for all times. 

As found above, the system's memory in this case is not sufficiently characterized by df. 
The mean coordination number z of the fractal subset, defined as the spatial average of the 
number of active neighbour sites, is not adequate either for the description the dynamics or the 
steady state. Indeed, the fractal set obtained by the iteration of the segment 1110011110011000 
has the same value of df and z = | as set Ci, but its dynamics and steady state are not the 
same due to the different fc-tuplet distribution. 

An alternative analytical description for the reaction on linear sets is provided by the 
theory of Markov chains. We have seen that the reactive dynamics of, say, set C%, can be 
fully determined by knowledge of the evolution of a single three-particle island. An island 
can never evolve into the empty state (SSS). There are 7 possible states ( 2 ) of such an island 
with at least one particle, namely (AAA), (SAA),(ASA),(AAS),(SSA),(SAS) and (ASS). Let 
us denote them by 1,2,3,4,5,6 and 7 respectively. The state vector of the system at time t is 
P(t) — (-F\(t), • • • ,Pi(t)) T , where P n (t) is the probability that the island be in the state n at 
time t. P(t) satisfies the stochastic evolution equation 

P(t + At)=W T P(t), (5) 



( 2 )The number of relevant states can be decreased by symmetry considerations, but we shall keep all seven 
states for the sake of clarity. 
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where the elements Wij of the transition probability matrix W are computed from the MC 
algorithm for the reaction by counting the number of paths leading from one state to another: 



/ 1/3 1/6 1/3 1/6 



W = 







2/3 1/6 1/6 

1 

2/3 1/6 1/6 

1 

1 

yooooooi/ 



(6) 



(we have set kn = 1 for simplicity). This matrix possesses four absorbing states and is therefore 
nonergodic The steady states II are the eigenvectors of the matrix W T corresponding to 
the degenerate eigenvalue Ai = 1. The algebraic multiplicity of Ai is 4, equal to the dimension 
of the subspace spanned by the four eigenstates IT with the components IT j = S^ r j, H2J = 

5^j,U.3 t j = 8qj and H4J = 5*j, where j — 1, .. .,7. The general form of the steady state IT 
can be expressed as a superposition of the states IT, i.e. n = J^i °i 3» with ^\ Cj = 1. The 
coefficients Cj depend on the initial conditions and can be calculated numerically by iteration 
of the evolution equation (JSJ. Thus, if one starts with a full island (state 1), the asymptotic 
mean number of particles on the island will be (N (00)) i S i. = 2- | + l- | + l- | + l- i = §, i.e. 
we recover the result obtained previously for the Cantor set C\. (JV(oo)) can be computed 
similarly for the set C2 and for any other Cantor-like sets. 

The other eigenvalues of W T describe the decay of {N(t)) to the steady state, i.e. (N(t)) — 
(iV(oo)) = Y^i=2 u i e ~ Xit w ith time-independent coefficients 1*2,3. They are A2 = § (twice) 
and A3 = 3. This time dependence is corroborated by eq. J2J), except that the arguments 
in the exponential functions differ from our result by a factor of ^. However, this artefact is 
a direct consequence of the choice for the time unit: in the model for the fc-tuplets and in 
the MC simulations, it was the time needed to scan the whole lattice, while here it has been 
implicitly assumed to be the time required to update a single site in a three-site lattice. 

Dynamics on fractals: diffusing reactants. - Reactive events in the above diffusionless 
systems can be viewed as particle jumps into already occupied sites. One can also allow for 
additional diffusion events, i.e. particles can jump into the "0" sites of the fractal sets with 
a probability rate kn- The "0" sites can then no longer be considered as inactive sites but 
rather as initially empty sites. A particle initially placed on a " 1" site of the fractal set can 
now diffuse into a "0" (empty) neighbour site. The latter becomes then a "1" site (occupied), 
while the original site becomes "0" (empty). On long time scales, this gives classical diffusion. 
A description in terms of islands is no longer suitable, since they can now interact with each 
other by means of diffusing particles at the boundaries of each island. 

To study the long-time dynamics, we have performed a series of MC simulations. Initially, 
particles are put on the " 1" sites of the lattice and then they start to diffuse freely over the 
entire lattice and react with each other. In this case we take reflecting boundary conditions, 
i.e. when particles arrive at sites or L + 1, they bounce off and go back to sites 1 and L, 
respectively. Again, the time scale is set by L . The reaction rate kn, taken to be equal to 
the diffusion rate ko, can be again absorbed in the time scale by setting r = kiit. 

For comparison with the diffusionless case, we have used as initial conditions the fractal 
sets C\ and C2 and have performed the same number of runs (2 • 10 3 ) for L = 4 5 . The time 
evolution of the concentration 9(t) has been monitored for sets C\ and C2. Figure [3 displays 
a double-logarithmic plot of numerical MC results for 8(t) as a function of r. The dynamical 
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Fig. 1 - Time evolution of the particle concentration in the absence of diffusion. The dots and crosses 
represent the MC results, while the continuous and dashed curves have been computed from eq. @. 

Fig. 2 - Time evolution of the particle concentration in the presence of diffusion. The slope of the 
continuous line is equal to — d//2, i.e. the theoretical value of the dynamical exponent (see text). 



exponent is given by the slopes of the point curves, which are represented by dots for C\ and 
crosses for C2; for both sets, the mean slope is very close to the value —df/2 represented by 
the continuous straight line. In contrast to the immobile reactant case, df appears to be the 
natural parameter associated to the long time decay of the concentration. 

The above analysis suggests a power law behaviour 9{t) = a t~ d f/ 2 governing the long time 
dynamics. This law exhibits two types of memory effects: a memory of the dimensionality 
df of the initial particle distribution entering via the exponent; and a more detailed memory 
of the initial distribution, entering via the amplitude factor a. It is only in the regime of 
extremely long times t 3> L 2 that the decay curves for sets C% and Ci fall into each other. 
This limit is to some extent trivial from the standpoint of many-particle dynamics, since 
only one particle remains on the lattice. In the long time regime, detailed information about 
lacunae in the initial particle distribution will also be kept before reaching the steady state in 
higher order quantities like the distribution function for the interparticle distance. 

These results agree with recent theoretical predictions. When k D = k Rl a closed analytic 
description in terms of empty fc-tuplets (also termed intervals) is possible |11I16| . This method 
has been used to show that, for fractal initial conditions with df < 2, (N(t)} behaves as t~ dl ' I 2 
at long times in the large L limit |lfilll7| : for df = 2, logarithmic corrections are necessary, 
while MF theory applies for df > 2. In our case df — log(3)/log(4) rs 0.792, so the exponent 
"remembers" the initial particle distribution. 

Conclusions. - In the present work the reaction process A + A — ► A with fractal initial 
conditions has been studied both for immobile and diffusing reactants. In the diffusionless case, 
it was shown that the fractal dimension df does not suffice to characterize the dynamics and 
the steady state. The number of surviving particles at the steady state depends on the details 
of the initial distribution. In the presence of diffusion, multistationarity is suppressed and the 
steady state becomes universal and MF-like. However, a long tail characteristic of anomalous 
dynamics subsists, implying that the decay is governed by a power law (N(t)) oc t~ d rather 
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than an exponential. In contrast to the immobile reactant case, the memory of the initial 
condition is carried by the characteristic dynamical exponent d rather than by the steady 
state, where d is fully determined by df. As expected, the memory of the initial condition is 
less detailed than in the immobile reactant case, due to the randomizing effect of diffusion. 

The above results may be relevant for a series of experimental situations involving systems 
other than particle aggregates. A few examples are heterogeneous catalysis, evaporation- 
deposition systems, porous media, percolation clusters and ferromagnetic systems. 

Our work can be generalized in many different ways. One can e.g. consider the case of 
random rather than deterministic fractal initial conditions. Other possibility is studying more 
complex reactive schemes like the reversible case A + A ^ A 3 , whose long-time dynamics 
still remains to be characterized in detail in the limit of immobile reactants. 
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